Magnetic properties of PdAsaOg: a dilute spin system with an unusually high Neel temperature 
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The crystal structure and magnetic ordering pattern of PdAs206 were investigated by neutron powder 
diffraction. While the magnetic structure of PdAs206 is identical to the one of its isostructural 3d-homoIogue 
NiAs206, its Neel temperature (140 K) is much higher than the one of NiAs205 (30 K). This is surprising in 
view of the long distance and indirect exchange path between the magnetic Pd^"*" ions. Density functional 

calculations jaeld insight into the electronic structure and the geometry of the exchange-bond network of 
both PdAs206 and NiAs206, and provide a semi-quantitative explanation of the large amplitude difference 
between their primary exchange interaction parameters. 



PACS numbers: 75.25.-j, 75.50.Ee, 61.05.ftn, 71.20.Ps 

I. INTRODUCTION 

The magnetic properties of transition metal compounds 
with 3d valence electrons have been one of the central 
research themes in solid-state physics for the past three 
decades. In view of the interplay between magnetism and 
high-temperature superconductivity, particular attention 
has been focused on oxides and arsenides. Recently, the 
electronic structure and ordering phenomena of transition 
metal compounds with 4d and 5d valence electrons (such as 
ruthenates and iridates) have also captured much attention. 
The electronic correlations in these materials are generally 
weaker than those of their 3d counterparts, while the spin- 
orbit coupling is stronger. The quantitative description of 
the influence of these parameters on the electronic phase 
behavior of d -electron compounds is an important topic of 
current research. Here we report a detailed investigation of 
the magnetic properties of PdAs206, a recently synthesized^ 
electrically insulating compound with a magnetic lattice of 
Pd^"*" ions in the electron configuration 4d^. We compare 
our results to the isostructural compound NiAs206, which 
is based on Ni^"*" ions with the same number of electrons in 
the 3d -shell. 

PdAs206 crystallizes in the PbSbjOe structure with Pd^+ 
and As^""" ions segregated into different layers (Fig. 1). 
The octahedral coordination of Pd^"*" in this structure is 
unusual, because divalent palladium shows a strong pref- 
erence for square-planar coordination, which is associated 
with a diamagnetic ground state. Only a few examples of 
sixfold-coordinated Pd^"*" compounds are known, including 
the ambient- and high-pressure polymorphs of PdFj as well 
fluoro-palladates of composition MPdF4 (M = Ca, Cd, Hg) 
and CsPd2F5.'^ These compounds are paramagnetic at high 
temperatures and tend to order antiferromagnetically upon 
cooling. In accord with this trend, magnetic susceptibility 
measurements on PdAs206 showed paramagnetic behavior 
at room temperature, and an antiferromagnetic phase tran- 
sition at the Neel temperature T^^ ~ 150 K.^ This behavior 
is qualitatively analogous to the one of the isostructural 
3d-homologues MnAs206, C0AS2O6, and NiAs206, which 



also show antiferromagnetic ordering with Tjv = 13, 20, 
and 30 K, respectively.'^ However, the much higher Neel 
temperature of PdAs206 is surprising, especially because 
the PdOe-octahedra do not share vertices, edges or faces. 
The exchange paths connecting neighboring Pd^"*" ions are 
therefore long and involve at least two bridging oxygen 
sites. 

In order to elucidate the microscopic origin of this surpris- 
ing behavior, we have used neutron diffraction to determine 
the magnetic structure of PdAs205, which turned out to be 
identical to that of the 3d-homologues (Section II). This 
implies similar networks of exchange bonds in both sets of 
compounds. We employed density-functional calculations to 
obtain insights into the electronic structure of PdAs206 and 
NiAs206, and specifically into the origin of the exchange 
paths (which turned out to be hopping via As dimers) and 
of the large amplitude difference of the primary exchange 
interaction parameters (Section III). A model based on these 
interactions yields an excellent description of the magnetic 
susceptibilities of both compounds. 



II. NEUTRON DIFFRACTION 

A. Experimental details 

A powder sample of PdAsjOg of weight ~ 3 g was pre- 
pared using the starting materials PdO (99.9 % metals basis, 
Alfa Aesar) and AS2O5 (99.9 %, Alfa Aesar) in the molecular 
ratio 1 : 1.1 as described earlier.^ The mixed powder was 
pressed into pellets and dried in evacuated silica tubes for 
12 h at 373 K. Then the evacuated silica tubes were heated 
up to 973 K with a rate of 100 K/h. The hygroscopic and 
air sensitive powder of PdAs205 was obtained after an an- 
nealing process of about 100 hours. Measurements of the 
magnetic susceptibility and specific heat were carried out in 
the temperature range between 5 and 300 K (Fig. 2). Both 
quantities show anomalies indicative of antiferromagnetic 
ordering of the Pd sublattice at 140 K, in agreement with 
prior work.^ 




Fig. 1: (Color online) Crystal structure ofylASaOg (A = Pd, Ni) 
showing the three-dimensional network of AOg and AsOg octahedra. 
The brown (dark grey) colored balls denote A atoms, and the cyan 
(light grey) colored balls represent As atoms. Small black colored 
balls at the corners of the octahedra are O atoms. 
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Fig. 2: Temperature dependence of the specific heat and the inverse 
magnetic susceptibility of PdAs206. 



In order to investigate the crystal and magnetic structure 
of PdAsjOs, a neutron diffraction experiment was carried 
out at the research reactor FRM-II in Garching. Neutron 
powder patterns were collected with the instrument SPODI 
at 5 K and 200 K in the 20 range 4°-160°. This instrument 
uses a germanium monochromator (reflection 551) selecting 
the neutron wavelength A = 1.5476 A. The refinements of 



the crystal and magnetic structure were carried out with the 
program FullProf."* We used the nuclear scattering lengths 
b(Pd) = 5.91 fm, fa(As) = 6.58 fm and b(0) = 5.805 fm.'^ 
The magnetic form factors of the magnetic ions were taken 
from Ref.6. 



B. Crystal structure of PdAs206 

The crystal structure of PdAsjOg was recently refined 
from x-ray powder diffraction data in the trigonal PbSbjOs- 
type structure (space group P31m, No. 162), where the Pd, 
As and O atoms are in the Wyckoff positions la (0,0,0), 
2d(i,|,|) and 6fc(x,0,z), respectively.^ The same space 
group was found earlier for the compounds MnAsjOg, 
C0AS2O6 and NiAsjOg containing 3d-metal ions.^ From our 
neutron powder diffraction data taken at the lower temper- 
atures 5 and 200 K (Fig. 3) the trigonal space group P31m 
was confirmed. For the Rietveld refinements we used data 
in the extended 20 range from 4° up to 146°. A total of 14 
parameters was refined: an overall scale factor, five profile 
function parameters, the zero point, two lattice constants, 
the positional parameters x and z of the oxygen atom as 
well as three isotropic thermal parameters. The powder 
sample contained small amounts of the binary oxide PdO, 
which crystallizes in the tetragonal space group P42/mmc.^ 
Therefore the overall scale factor of PdO was additionally 
allowed to vary during the refinements. 

In Table I the results of the refinements are compared 
with those of the x-ray study carried out earlier at room 
temperature.^ Here it can be seen that the positional param- 
eters of the oxygen atoms determined at 5 and 200 K are 
in good agreement, indicating that the structural changes 
between the magnetically ordered and the paramagnetic 
states are weak. Only a slight reduction of 0.0036 A (about 
6 cr) could be observed for the Pd-O-bond length in the 
PdOg-octahedra. In contrast, the distances between the As 
and O -atoms are practically unchanged (Table I). The value 
d(As-O) = 1.8281(6) A (at 200 K) found for PdASjOg is 
in very good agreement with the values of other arsenates 
containing 3d-metal ions: d(As-0) = 1.827(4) A (NiAsjOg), 
d(As-O) = 1.830(3) A, (C0AS2O6), d(As-O) = 1.826(2) A 
(MnAs206), and d(As-O) = 1.826(1) A (CdASjOg).^-*^ All 
of these values are in agreement with d(As-0) = 1.82 A 
calculated for an AsOg-octahedron given by Shannon.' In 
Table I it can be seen that the structural parameters ob- 
tained at 200 and 300 K show relatively large discrepancies, 
despite the fact that both data sets were collected in the 
paramagnetic phase. This can be ascribed to the larger scat- 
tering power of the O -atoms in neutron diffraction, with the 
consequence that the 0-positions can be determined more 
reliably. Furthermore, the shortest oxygen contact d(O-O) 
= 2.308(3) A was found to be implausibly short in the x-ray 
study. ^ From our neutron diffraction study we found the 
larger values d(O-O) = 2.3757(12) A at 5 K and d(O-O) = 
2.3726(12) A at 200 K, respectively The value d(O-O) = 
2.410(3) A found for NiAs206 is shghtly larger,^ while the 
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Fig. 3: (Color online) Neutron powder diffraction data of PdAs206 collected at 5 and 200 K. The crystal structure was refined in the 
hexagonal space group P31m. The calculated patterns (red lines) are compared to the observed ones (black circles). In the lower part of 
each diagram the difference pattern (blue) as well as the positions of the nuclear reflections of PdAs205 and the impurity phase PdO 
(black bars) are shown. For the powder pattern collected at 5 K, the positions of the magnetic reflections of PdAs205 are also shown. The 
inset shows the magnetic Bragg reflections of PdAs205, obtained from the difference between the diffraction patterns at 5 and 200 K. 



As-O-bond lengths d(As-O) = 1.827(4) A (NiAs206) and 
d(As-O) = 1.8281(6) A (PdAs206) are practically the same 
in both compounds. The cell volume of the Ni compound 
(V = 86.97(3) A^) is much smaller than the one of the Pd 
compound (V = 93.698(2) A^). This is due to fact that 
the ionic radius of Pd^"*" is larger than that of Ni^"*". In or- 
der to keep the As-O-bond lengths almost constant in the 
AsOe-octahedra, the bond angle Z(O-As-O) increases from 
169.91(3)° in PdASaOg) to 173.23(12)° in NiASjOfi. 



C. Magnetic structure of PdAs205 

The neutron powder data recorded at 5 K show additional 
Bragg reflections that can be ascribed to the antiferromag- 
netic order of the Pd sublattice. The two prominent ones at 
26 = 9.5° and 26 = 23.5° can be indexed as (0, 0, |)m and 
(0, 1, |)j^, respectively. This suggests that the magnetic cell 
has a doubled c-axis with a propagation vector k = (0, 0, 
i). All further magnetic reflections were assigned indices 
according to {hk£)^ = QiM)^^ ± k, where M and N desig- 
nate the magnetic and nuclear reflections. The same type 
of magnetic ordering was observed earlier for the isotypic 
divalent transition metal arsenates C0AS2O5 and NlAs206.^ 
The presence of the strong magnetic reflection (0, 0, ^)jvf 
indicates that the magnetic moments of the Pd-atoms are 
aligned ferromagnetically within the hexagonal ab -plane. 



Due to antiferromagnetic exchange interactions between the 
palladium moments the ferromagnetic layers form the se- 
quence -I 1- - along the c-axis (Fig. 4). With this model the 

magnetic structure of PdAs206 could be successfully refined 
using the magnetic reflections observed in the 20 -range 
up to 50°. It has to be noted that the moment direction 
within the hexagonal plane cannot be determined from the 
refinements. The wave vector k = (0, 0, i) keeps the full 
symmetry of the group G^ = G according to magnetic group 
theory, and it defines the magnetic translation lattice.^" The 
existence of three magnetic domains in the hexagonal ba- 
sis plane prohibits an unambiguous determination of the 
moment direction. 

In order to improve the refinement of the magnetic struc- 
ture, we used the purely magnetic intensities obtained 
from the difference between the data sets collected at 5 
and 200 K (inset of Fig. 3). Since the magnetic form 
factor of Pd^""" is not available, we first used the one of 
the Pd"*" ion,^ but Table II shows that the calculated in- 
tensities decrease much more strongly with increasing 26 
than the observed ones. A considerably better fit was ob- 
tained with the form factor of Ni^"*", which also shows a 

-configuration. The relatively large residual Rj^ = 0.141 
[defined as = Y^WOobsl - \Icaic\\/\Uobs\] reflects the fact 
that the low intensities of the magnetic reflections at high 
20 -values could not be measured with good accuracy, in 
combination with systematic errors arising from the differ- 
ence between the form factors of Pd^+ and Ni^"*". 
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Table I: Results of Rietveld refinements of the neutron powder 
diffraction data (A = 1.5476 A) for the nuclear structure of 
PdAs206 at 5 and 200 K. The lattice constants, positional and 
isotropic thermal parameters as well as the bond distances and an- 
gles within the AsOg- and PdOg-units are compared with the values 
obtained earlier at room temperature from x-ray powder diffrac- 
tion data (A = 0.7093 A).^ The residual Rjv of the refinement of 
the crystal structure is defined as R„ = Z IK^otJ - l^'ccfcll/KfoiJ- 
The residual for the room temperature structure (marked by *) 
was calculated with intensities rather than structure factors.^ 
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Table II: Observed and calculated intensities of the magnetic re- 
flections of PdAs206 as obtained from Rietveld refinements using 
the magnetic form factors of Ni^""" and Pd^.^ The residual is 
defined as = S llUobJ - \haic\\/\Qobs\- 



The sublattice magnetization resulting from the refine- 
ment is /i = 1.92(4)|U5 per Pd site, similar to the value 
/i = 2.11(3)/^B reported for NiAsjOg (Ref.3). While the 
ordered moment is consistent with the spin-only moment 
expected for a configuration, a fit to the magnetic sus- 
ceptibility for T > Tj^ yields a g-factor larger than 2, which 



is indicative of an orbital contribution to the Pd moment 
(see Section III.C below). The difference may in part be due 
to zero-point fluctuations of the magnetic moment, which 
reduce the ordered moment of the spin-1 system in the 
binary oxide NiO by ~ 8 The zero-point reduction is 
possibly larger in PdAsjOg because of the low-dimensional 
exchange-bond network (see Section III.C). While these con- 
siderations suggest a small but nonzero orbital contribution 
to the Pd moment, measurements of the g-factor anisotropy 
by single-crystal neutron diffraction and/or electron spin 
resonance will be required to quantify this contribution. 




Fig. 4: (Color online) Magnetic structure of PdAs206. Shown 
are the isolated PdOg-octahedra between the AsOs-layers. The 
magnetic moments of the palladium atoms are aligned ferromag- 
netically within the hexagonal ab-plane. Along the c-direction 
the moments are coupled in an antiparallel fashion, forming the 
sequence H h -. 



III. DENSITY FUNCTIONAL CALCULATIONS 
A. LDA band structure 

Figure 5 shows the electronic band structure and den- 
sity of states in the paramagnetic local density-functional 
approximation (LDA)^^ for NLAsjOe in solid lines and for 
PdAsjOe in dashed lines. The self-consistent calculations 
were performed with the linear muffin-tin orbital (LMTO) 
method^'^ using 8x8x8 k-points in the Brillouin zone. 

The Fermi level falls in the middle of the two narrow 
transition-metal (^4) -bands. These are split above the 
three narrow tjg-bands, because the pdcr-hopping is ~ ^/3 
times stronger than the pdTi -hopping to the O p orbitals at 
the corners of the AOe-octahedron. This is in accord with 
the labeling Similarly, in accord with the labeling 
and As^+, the oxygen 2p-like bands are below and the As s- 
and As p-like bands are above the Fermi level. However, in 
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Fig. 5: (Color online) LDA band structure (left panel) and density 
of states (right panel) for NLAS2O6 (solid lines) and for PdAs206 
(dashed lines) . The two band structures have their Fermi levels 
(0 eV) lined up. The bands are plotted along the high-symmetry 
directions of the hexagonal Brillouin zone and the densities of 
states have their dominating characters labeled. 



terms of atomically localized orbitals (LMTOs) rather than 
Wannier orbitals, the bands denoted as As s in the figure 
have ~40% anti-bonding O p character, as well as some A s 
and Os characters. Correspondingly, around — 12 eV (below 
the frame of the figure) there are two bands with As s and 
O sp bonding characters in about equal amounts. 

With the bands lined up at the Fermi level, the and 
p bands lie lower in the Pd than in the Ni compound. This 
is because 4d orbitals have a larger extent and therefore 
larger hopping integrals to O p than do 3d orbitals. For the 
same reason, the band is about 1.4 times wider in the Pd 
than in the Ni compound. 

Inclusion of the Coulomb interaction beyond the LDA 
splits the Cg bands and leads to insulating solutions, as 
we have checked through LDA+U^'^ calculations. For the 
present purpose of calculating and understanding the mag- 
netic properties, it is more convenient to start from the 
localized description and treat the hopping, t, to order 



B. Pdfig Wannier orbitals, low-energy tight-binding 
Hamiltonian and magnetic interactions 

We therefore construct a low-energy Hubbard Hamil- 
tonian. Since the LDA band is narrow and well sepa- 
rated from all other bands, we can limit the one-electron 
Hilbert space to that of the two A-centered Wannier or- 
bitals, 3z^ — and — y^, describing this band. When 
using the iVMTO downfolding method, we thus kept the 
A Eg degrees of freedom as active, and downfolded the 
rest. The Wannier orbitals are finally obtained by sym- 
metrically orthonormalizing the downfolded, iVth-order 
muffin-tin orbitals (JVMTOs).^^ In the representation, H = 
S '^im,jm'(.^im^jm' + h.c), of the Corresponding one-electron 
part of the Hamiltonian, is the hopping integral from 

orbital m on site i to orbital m' on site j. 

Figure 6 shows the — and 3z^ — Wannier orbitals 
for PdAsjOs as x (r) = ±const. surfaces with positive lobes 
colored red and negative blue. The z-axis points along the 
Pd-0 bond in the bc-plane and the x- and y-axes point to the 
other oxygens in the Pd-centered octahedron. In the figure, 
the latter has been given a blue, transparent skin while all 
As-centered octahedra have been given a red, transparent 
skin. Pd atoms are green. As atoms yellow, and oxygens 
violet. Such a Pd-centered Wannier orbital has — y^ 
or 3z^ — character locally, and strong pda antibonding 
character on the neighboring oxygens. The unusual feature 
here is that the bacfc-lobes of the strongest O p characters 
(the four p^ and py tails of the x^ — y^ Wannier orbital and 
the two tails of the 3z^ — Wannier orbital) bond to 
the sp characters on the closest pair of As atoms and, from 
there, antibond to the closest oxygen, which now belongs 
to a neighboring PdOe octahedron. As an example: The 
red lobe of the x^ — y^ orbital sticking up and out towards 
the reader, antibonds with the blue lobe of near p^ orbital, 
whose red back-lobe merges together (bonds) with the red, 
two-center As sp bond, giving rise to a "red sausage". The 
latter finally antibonds with the O p orbital which points 
upwards towards a near Pd belonging to the upper Pd-sheet 
(not shown in the figure). Similarly for the 3z^ — orbital 
on the right-hand side of the figure: Its red lobe, sticking 
down and out, antibonds with the blue lobe of the near p^ 
orbital, which merges with the As two-center bond into a red 
sausage. The latter finally induces strong -like character 
which points down and out, towards a Pd atom in the lower 
Pd sheet. 

These hopping paths, shown schematically in Fig. 7, 
are to 3rd-nearest Pd neighbors, but only to six out of 

the twelve, namely to those at ± (0, 1,1), ± ^, -l] , 
and i ("Y' ~2' ''^} ■ integrals for hopping to the 

3rd-nearest neighbors at ±(0,1,-1), ±(^,-,1), and 

± ^ — ^, 1 j , which have no bridging oxygen and As pairs, 
are negligible. The calculated hopping integrals exceeding 
10 meV are given in Table III. We see that those to 3rd- 
nearest neighbors dominate those to 2nd- and Ist-nearest 
neighbors. 
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Fig. 6: (Color online) The — (left panel) and 3z^ — (right panel) Wannier orbitals which span the LDA Pd bands of PdAs206. 
Plotted are orbital shapes (constant-amplitude surfaces) with lobes of opposite signs colored red (dark) and blue (light) . For clarity, the 
central PdOg octahedron is given a blue skin and all AsOg octahedra a red skin. Pd atoms are green, As atoms yellow, and oxygens violet. 




Fig. 7: (Color online) Paths of the dominant hoppings. Dashed 
green: between — orbitals (see left-hand side of Fig. 6) on 
3rd-nearest A neighbors. Solid red: between 3z^ — orbitals (see 
right-hand side of Fig. 6) on 3rd-nearest A neighbors. The colors 
of the atoms are as in Fig. 1. 



The Ni 3d [e^j Wannier orbitals in NiAs206 are similar 

to the Pd 4d {e^ Wannier orbitals in PdAsjOe, except that 
they are more localized. This is consistent with the 1.4 times 
smaller 3d [e^-hand width. Concomitantly, the hopping 
integrals for NiAsjOe listed in parentheses in Table III are 
about 1.4 times larger than those for PdAs205. 

The magnetic exchange interaction, J, can be expressed 
in general as a sum of antiferromagnetic and ferromagnetic 



contributions. In the limit of large Coulomb correlation, 
typically valid for late transition metal elements like Pd 
or Ni, the antiferromagnetic contribution is related to the 
hopping integral t by the second-order perturbation relation 
~ Y; ' where U is the effective on-site Coulomb repulsion 
defined for the Wannier orbitals. Considering the hopping 
integrals in Table III, the contribution from the term in 
the magnetic exchange gives rise to a factor two difference 
between the Pd and Ni compounds. The estimate of U, 
unlike that of the hopping integral t, is a rather delicate 
issue. In principle, one should compute U for the Wannier 
orbitals shown in Figure 7. However lacking a prescription 
to do so, we computed the U values corresponding to Ni- 
d and Pd-d partial waves, which were truncated outside 
the atomic spheres defined around an Ni or Pd sites. The 
calculations were carried out within the framework of the 
constrained DFT scheme. The U values calculated in 
this way were Up^ = 5.7 eV and t/jv; = 8.4 eV But since 
the Wannier orbitals are far more delocalized than the 
truncated partial waves, and more so for Pd than for Ni, 
these U values should be significantly reduced, and more 
so for Pd than for Ni. This could conceivably lead to the 
factor Uj^JUp^ ~ 2.4 needed to bring our theoretical t^/U 
estimate for the Neel temperatures of the two compounds 
into agreement with the measured ratio of 4.7. 



C. Spin model and susceptibility 

Taking into account only the dominant hopping integrals 
provided by the iVMTO-downfolding calculation, a spin 
model can be defined in terms of the six pairs of 3rd-nearest- 
neighbor magnetic interactions, all of size J3, obtained by 
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Table III: Hopping integrals between Wannier orbitals centered on Pd (or Ni). Hoppings are in meV and listed are those exceeding 10 
meV Only hoppings between 2nd-nearest neighbors, 4.82 (4.76) A, and between 3rd-nearest neighbors, = 6.71 ( 6.51) A, are 
significant. 



summing over the squares of the hopping integrals be- 
tween the 3rd-nearest neighbors. Based on this model, we 
have calculated the magnetic susceptibihty by considering 
the following spin-1 Hamiltonian: 



m— 1 m— 1 m— 1 

^ ~ '^3 ^ ^ ^ (^i,j,i<^L,)+l,ic+l + ^i,;,)c'^i+l.j.fc+l 
k=0 i=0 i=0 

This model was solved by the quantum Monte Carlo 
method^' on a 10 x 10 x 10 lattice. The primary inter- 
action J3, and the effective Lande g-factors were obtained 
by fitting to the experimental susceptibility. As shown in 
Fig. 8, the calculated and measured susceptibilities agree 
very well. The optimal values of the g -factor and the ex- 
change parameter J3, were found to be respectively 2.38 
and 62 K for the Pd compound, and 2.48 and 13 K for the 
Ni compound. We find that the g-factors are larger than the 
spin-only value of 2, in accord with the discussion in Section 
II. C above. The fit of the susceptibilities yields a value of 5 
for the ratio of the predominant magnetic exchange param- 
eters of PdAsjOe and NLAsjOe, in good agreement with the 
ratio of the magnetic transition temperatures. 



IV CONCLUSIONS 



Using a combination of susceptibility and neutron diffrac- 
tion measurements, we have developed a comprehensive 
experimental description of the magnetic properties of the 
newly synthesized antiferromagnet PdAsjOe. Density func- 
tional theory has provided a detailed understanding of the 
magnetic bond network of this compound, as well as a semi- 
quantitative explanation of the large enhancement of the 
magnitude of its primary exchange interaction parameters 
compared to its 3d homologue NiAsjOg. This approach may 
prove useful for research on other Pd compounds including 
the recently discovered^" ferromagnet PdS207, and for com- 
parative studies of materials with 4d valence electrons and 
their 3d -electron counterparts in general. 




T(K) 



T(K) 



Fig. 8: (Color online) Temperature dependence of the magnetic 
susceptibility (multiplied by temperature) for NiAs206 (left panel) 
and PdAs205 (right panel). The circles correspond to experimental 
data (taken from Ref.3 for NiAs206), and the sohd lines correspond 
to calculated data based on the model described in the text. 
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